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^ ; Abstract 

We consider the flow of an electrically conducting fluid between differentially rotating cylinders, 
in the presence of an externally imposed toroidal field Bo(fi/r) e^. It is known that the classical, 

>. 

axisymmetric magnetorotational instability does not exist for such a purely toroidal imposed field. 

m ' 

We show here that a non-axisymmetric magnetorotational instability does exist, having properties 
very similar to the axisymmetric magnetorotational instability in the presence of an axial field. 
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The magnetorotational instability (MRI) is a mechanism whereby a hydrodynamically 
stable differential rotation flow may be destabilized by the addition of a magnetic field. An 
obvious question then is whether different orientations of the field yield different types of 
instability, or the same, or none at all. The original view was that the axial component of 
the field is the only important one, with any azimuthal component playing no essential role, 
and incapable of producing any instabilities on its own [l|, Q] . This view was altered by the 
discovery that a mixed axial and azimuthal field yields instabilities quite different in many 
ways from those found with a purely axial field |3j, IJ]. In this letter we show that even a 
purely azimuthal field yields an MRI, and compare its properties with the previously known 
types. 

While its most important application is to astrophysical accretion disks [2J, the MRI was 
originally discovered in the much simpler Taylor- Couette problem, consisting of the flow 
between differentially rotating cylinders [l| • Because of its relative simplicity, this geometry 
has proven particularly amenable both to theoretical analyses of the MRI m, 3] , as well as 
to its recent experimental realization j^, We too shall examine the MRI in this context. 

Consider therefore an electrically conducting fluid confined between two concentric cylin- 
ders of radii rj and r G , rotating at rates fli and Q , chosen to satisfy fi r o > ^i r i- That is, 
the angular momentum increases outward, so by the familiar Rayleigh criterion the flow is 
hydrodynamically stable, with the angular velocity given by 

Q{r) = A + B/r 2 , (1) 

where 

' o i o i 

In this work we will fix r a = 2rj and Q Q = fij/2, so A and B simplify to |fij and f^rf , 
respectively. The essence of the MRI then is to ask whether the addition of a magnetic field 
can destabilize this flow. 

If the imposed field is purely axial, B = B e z , the profile (1) can be destabilized, 
provided the rotation rates are sufficiently great, and the field strength B is neither too 
weak nor too strong; . Specifically, the magnetic Reynolds number Rm = Qirf/rj, where 
r] is the magnetic diffusivity, must exceed O(10), and the Lundquist number S = B ri/r]^/pJI, 
where p is the density and \i the permeability, must be around 3 — 10. 
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In contrast, if the imposed field is mixed axial and azimuthal, Bo = B$ e z + fiB^irijr) e^, 
where (3 is around 1 — 10, then theprofile (1) can again be destabilized, but at very different 
rotation rates and field strengths |3j, |4j. The relevant parameter measuring the rotation rates 
is now not the magnetic Reynolds number, but rather the hydrodynamic Reynolds number 
Re = Qirf/u, where v is the viscosity. Similarly, the parameter measuring the field strength 
is not the Lundquist number, but instead the Hartmann number Ha = B Q rij J~pin\v. The 
MRI sets in when Re > O(10 3 ), and Ha » O(10). 

To compare these results, we note that the two sets of parameters are related by 
Re = RmPm -1 and Ha = SPm _1//2 , where Pm = vji] is the magnetic Prandtl number, 
a material property of the fluid. Typical values for liquid metals are O(10 -6 ). Translating 
the results for the purely axial field, we thus obtain Re ^ O(10 7 ) and Ha ~ O(10 4 ), both 
several orders of magnitude greater than for the mixed field. It is perhaps not surprising 
then that the MRI has been obtained experimentally for the mixed field 

as 

, but not (yet) 

for the purely axial field [5(. 

As different as they are, one feature these two types of MRI have in common is that 
they are both axisymmetric. Non-axisymmetric modes have also been explored, for both 
the purely axial p| as well as the mixed fields |4j]. For a purely axial field the relevant 
parameters are still Rm and S, but the critical values Rm c are somewhat larger than for 
the axisymmetric modes, indicating that the axisymmetric MRI is the most unstable mode. 
For mixed fields, one finds — perhaps somewhat surprisingly — that adding an azimuthal 
component now has minimal effect, certainly far less than the reduction by four orders of 
magnitude found for the axisymmetric modes. Evidently the relevant parameters continue 
to be Rm and S, rather than Re and Ha. 

What we wish to show in this letter then is that for these non-axisymmetric modes, one 
can impose a purely azimuthal field, B (ri/r) e^, and still obtain an MRI, having all the 
characteristics of the previous non-axisymmetric types of MRI 0, Q] . To this end, we solve 
the linear stability equations 

db 

Rm— = V 2 b + V x (u x B ) +RmV x (U x b), (3) 

Re = -Vp + V 2 u + Ha 2 (V x b) x B + Re (U x V x u + u x V x U ), (4) 
at 

where Uo = rVt{r) is the profile (1) whose stability we are exploring, and B = B (ri/r) 
is the imposed azimuthal field. Length has been scaled by r i; time by , U and u as 



r^i, Bo as Bq, and b as RulBo- 

Taking the t, z and <fi dependence to be of the form exp(crt + ikz + imtfi), and using also 
V • b = to eliminate b z , the r and <fi components of (3) become 

Kmab r = V 2 6 r — r~ 2 b r — 1imr~ 2 b$ + imr~ 2 u r — KmimQb r , (5) 

Rm ab^ = V 2 b ( / > — r~ 2 b^ + 2imr~ 2 b r + imr~ 2 u^ + 2r~ 2 u r — RmimQb^ + RmQ'rb r , (6) 

where primes denote d/dr. The components of (4) have a similar structure, but we will 
not need to refer to them in the subsequent discussion, and hence do not list them. The 
boundary conditions associated with (3) and (4) are 

b r = b' (j) + r^btp = u r = = u z = (7) 

at r = ?~j and r , corresponding to perfectly conducting, no-slip walls. The resulting one- 
dimensional linear eigenvalue problem is solved by finite differencing in r, as in [7J]. 

Figure 1 shows stability curves for m = 1, the only mode that appears to become unstable. 
We see how an MRI exists that is remarkably similar to some of the results described above. 
In particular, as Pm — > 0, the relevant parameters are clearly once again Rm and S, with 
the MRI arising if Rm > 80, and S ~ 40 yielding the lowest value of Rm c . The specific 
numbers are roughly an order of magnitude greater than for the axisymmetric MRI in the 
axial field, but the basic scalings, and even the detailed shape of the instability curves, are 
identical. 

Figure 2 shows the real and imaginary parts of o in the unstable regime. Remembering 
that time has been scaled by fiT" , we see that we obtain growth rates as large as 0.05fij. So 
again, while the particular number 0.05 is about an order of magnitude smaller than for the 
axisymmetric MRI in the axial field, this non-axisymmetric MRI is clearly also growing on 
the basic rotational timescale. 

To understand why this non-axisymmetric MRI exists even for a purely toroidal field B , 
for which it is known that the axisymmetric MRI fails l|, |2[], we need to consider the details 
of (5) and (6). In particular, note that for m = 0, b r completely decouples from everything 
else, and inevitably decays away. Without b r though, the MRI cannot proceed, as it relies 
on the term R,mQ'rb r in (6). In contrast, for m = 1, b r is coupled both to b^, coming from 
V 2 b, and to u r , from V x (u x B ). And once b r is coupled to the rest of the problem, 
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FIG. 1: The critical magnetic Reynolds number for the onset of the MRI, as a function of the 
Lundquist number, for the different values of Pm indicated, m = 1, r /ri = 2, £l /Qi = 0.5. 
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FIG. 2: The real (top) and imaginary (bottom) parts of a, as functions of S, with Rm fixed at 200. 
The real part corresponds to the growth rate, the imaginary part to the azimuthal drift rate. 

the term RmQ'rb r then allows the MRI to develop. Figure 3 presents an example of these 
solutions, indicating how all three components of both u and b are indeed present. 

There is however one aspect that is not entirely clear from this analysis, namely why 
this non-axisymmetric MRI actually requires the term RmQ'rb r at all. In particular, the 
axisymmetric MRI with a mixed field does not rely on it, but instead on the term 2r~ 2 u r , 
coming from V x (u x B ) rather than RmV x (Uo x b) [3j. The non-axisymmetric MRI 
is evidently more like the axisymmetric MRI with a purely axial field, which also requires 
the term RmQ'rb r , because in that case the term 2r~ 2 u r is absent. 
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FIG. 3: The marginally stable solution at Pm = 0.1, S = 47.4 and Rm = 93 (see Fig. 1). From left 
to right the r, 4> and z components of u (top) and b (bottom). The real parts are solid, imaginary 
parts dashed. Note how both u and b have the z components largest. The azimuthal wavenumber 
is k = 1.88. 

To summarize then, we have shown that even if the externally imposed magnetic field is 
purely toroidal, one still obtains an MRI, simply non-axisymmetric rather than axisymmet- 
ric. In other respects though this new instability is remarkably like the classical axisymmetric 
MRI with a purely axial field, indeed far more like it than the axisymmetric MRI with a 
mixed field, which yielded fundamentally different scalings. In contrast, the parameters here 
continue to be Rm and S, just as in the classical MRI. 

Note though that attempting to obtain this non-axisymmetric MRI in a laboratory ex- 
periment would be even more difficult than attempting to obtain the classical axisymmetric 
MRI in an axial field. First, the re quir ed rotation rates would be even greater, Re > O(10 8 ), 
with all the difficulties that entails [10]. Even more daunting, imposing an azimuthal field of 
the required strength, Ha ~ O(10 5 ), would require a current along the central axis in excess 
of 10 6 A, surely far beyond any feasible experiment. 
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This new non-axisymmetric MRI could have astrophysical applications though, since 
many astrophysical objects do have predominantly azimuthal fields, in which case the re- 
sults presented here suggest that this non-axisymmetric MRI could be preferred over the 
axisymmetric MRI. 

Despite the similarities in their fundamental scalings, the fact that the classical MRI is 
axisymmetric, whereas this new MRI is non-axisymmetric, means the nonlinear equilibra- 
tion, and hence the associated angular momentum transport, are potentially quite different, 
which could again have astrophysical implications. Work on this is currently in progress. 

Finally, one might ask how the results presented here change if one allows for a more 
general toroidal field profile, B = (cir" 1 + C2r) e^, where the term c^t corresponds to an 
electric current flowing through the fluid itself, and not just along the central axis. Eq. 
(4) then contains an additional term Ha 2 (V x Bo) x b, which opens up the possibility of 
instabilities driven entirely by this current V x B , without any rotation necessarily present 
at all HI]. Understanding how these current driven instabilities (also m = 1) interact with 
the magnetically catalyzed but ultimately rotationally driven MRI presented here is also in 
progress. 

This work was supported by the German Leibniz Gemeinschaft, under program SAW. 
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